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ABSTRACT 

We present ray tracing simulations combined with sets of large N-body simulations. 
Experiments were performed to explore, for the first time, statistical properties of 
fluctuations in angular separations of nearby light ray pairs (the so-called lensing 
excursion angle) induced by weak lensing by large-scale structures. We found that the 
probability distribution function (PDF) of the lensing excursion angles is not simply 
Gaussian but has an exponential tail. It is, however, found that the tail, or more 
generally non-Gaussian nature in the PDF has no significant impact on the weak 
lensing of the CMB. Moreover, we found that the variance in the lensing excursion 
angles predicted by the power spectrum approach is in good agreement with our 
numerical results. These results demonstrate a validity of using the power spectrum 
approach to compute lensing effects on the CMB. 
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1 INTRODUCTION 

Weak lensing effects on the temperature anisotropy in the 
cosmic microwave background (CMB) have been recognized 
as a powerful probe of cosmology (see Mellier 1999; Bartel- 
mann & Schneider 2001 for reviews). Dark matter distribu- 
tion along the line of sight between the last scattering surface 
(LSS) and us deflects the light ray trajectories and induces 
distortions in the pattern of the CMB anisotropies. Since the 
gravitational lensing is directly sensitive to the matter dis- 
tribution up to the LSS, lensing signatures imprinted on the 
CMB may provide important information about the mat- 
ter distribution on large scales and at high redshifts. In this 
point of view, various methods have been proposed (Seljak 
1996; Bernardeau 1997; 1998; Metcalf & Silk 1997; Zaldar- 
riaga & Seljak 1998; 1999; Zaldarriaga 2000; Suginohara, 
Suginohara & Spergel 1998; Van Waerbeke, Bernardeau & 
Benabed 2000; Takada, Komatsu & Futamase 2000; Takada 
& Futamase 2001). These lensing signatures are generally 
small but are measurable with two planned satellite mis- 
sions, MAPQ and Planckf], and can help to break some of 
the parameter degeneracies in the CMB (Bernardeau 1997; 
Metcalf & Silk 1998; Takada & Futamase 2000). The change 
in a separation angle of two nearby light rays caused by weak 
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lensing (lensing excursion angle) plays a key role in studying 
weak lensing of the CMB, especially weak lensing effects on 
two-point statistics of the CMB. 

The analytical prediction of the weak lensing effect on 
the CMB power spectrum in modern cosmological models 
was first developed by Seljak (1994) based on linear per- 
turbation theory (the so-called power spectrum approach). 
While, in his subsequent paper, Seljak (1996) examined ef- 
fects of the nonlinearity in the density through analytic fit- 
ting formulae (Peacock & Dodds, 1996), and pointed out 
that the nonlinearity is very important on sub-degree scales. 
Furthermore statistics of the lensing excursion angles due to 
weak lensing are frequently assumed to be Gaussian without 
a rigorous basis (e.g., Seljak 1996). Numerical simulations 
are, therefore, needed for testing the validity and limitation 
of the semi-analytic approach. 

The purpose of this paper is to examine the statistical 
properties of the lensing excursion angles using ray-tracing 
simulations combined with large A^-body simulations for the 
first time. Since the simulations were originally constructed 
for the cosmic shear statistics (Van Waerbeke et al. 2001; 
Hamana et al. 2001b) the maximum redshift is taken by z ~ 
3, which is about a half way to the LSS. This choice may not 
seem to be enough for studying the weak lensing effects on 
the CMB. However, most of the contributions to the CMB 
lensing come from structures at z < 3. To take an example 
shown by Suginohara et al. (1998), in the COBE normalized 
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Table 1. Cosmological parameters. 
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SCDM model, the structures within z < 3 contribute 85% 
to the variance of the lensing angular excursion from the 
LSS (see Figure lb of their paper). The remaining 15% will 
not significantly alter the basic properties of weak lensing of 
the CMB, and therefore we can, at least, study its essential 
features. It should be, however, noticed that this 15% is 
surely the most interesting for probing the early stage of 
the large-scale structure formation with weak lensing of the 
CMB. 

The outline of this paper is as follows. Our models and 
method of the ray-tracing simulation are described in §2. 
In §3, we show the numerical resolution of our simulation 
and discuss its limitation. Using results of the ray-tracing 
simulations, we examine statistical properties of the lensing 
excursion angles in §4. We conclude in §5. 

Throughout this paper, we work in the comoving coordi- 
nates system. The cosmological parameters are denoted with 
usual notations; Hubble constant, Ho = 100/ikm/s/Mpc; 
the density parameter, Q m ; the cosmological constant, Qa; 
the variance of the density fluctuation in a sphere of radius 
8/i _1 Mpc, erf. 



2 NUMERICAL METHODS AND MODELS 

In this sec tion, we describe methods of the T V-b ody simula- 



tion (§2.1) and the ray-tracing simulation (§2.2), and sum- 



marize our models. 

2.1 iV-body simulation and the tiling technique 

The JV-body simulation data set were generated with a PM 
code (see Hamana et al. 2001b for a detail description). Each 
iV-body experiment involves 256 2 x 512 particles in a peri- 
odic rectangular box of size (L, L,2L). The mesh used to 
compute the forces was 256 2 x 512. The initial conditions 
are generated adopting the transfer function of Bond & Ef- 
stathiou (1984) with the shape parameter F = fi m /i. We 
adopt three cosmological models; two flat models with and 
without the cosmological constant and one open model. The 
amplitude of the power spectrum is normalized by the clus- 
ter abundance (Eke, Cole & Frenk 1996; Kitayama & Suto 
1997). In Table |l|, cosmological parameters in each model 
are summarized. 

A light-cone of the particles was extracted from each 
simulation during the run as explained in Hamana, Colombi 
& Suto (2001a). Our aim was that the light-cone covers 
a large redshift range, < z ^ 3, and a field of view of 
5x5 square degrees with a wide angular dynamic range. To 
do that, we adopted the tiling technique first proposed by 
White & Hu (2000). We performed AW = 11, 12 and 13 
independent simulations for SCDM, OCDM and ACDM, re- 
spectively, covering adjacent redshift intervals [«f" n , zf iax ], 
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Figure 1. Tiling configuration of A^-body simulation boxes for 
SCDM model. Dashed lines (cone) show the comoving angular 
diameter distance of ±2.5 arcmin. Note that the scale of the ver- 
tical axis is exceedingly enlarged for clarity, the long side of the 
rectangular boxes is parallel to the horizontal axis. 



i = l,... , A/box- The size of each simulation is chosen such 
that the portion of the light-cone in [z™ m , zf 13 *] (aligned 
with the third axis) fits the box-size as shown in Figure [ij 
This way, angular resolution is approximately conserved as 
a function of redshift, except close to the observer (to be 
discussed in §^|). Finally, in order to keep large scale modes 
that may contribute to the weak lensing, we impose the sup- 
plementary constraint L > 80/i~ 1 Mpc. As a resulF L follows 
the sequence with redshift summarized in Table 



2.2 Ray-tracing simulation 

Let us first summarize basic equations of the multiple lens- 
plane algorithm which are directly relevant to this paper 
(see e.g. Schneider, Ehlers & Falco 1992 and Jain, Seljak 
& White 2000 for details). We shall denote the comoving 
Cartesian coordinate system as (xi,X2,y) with y being the 
third axis and the origin being at the observer point. Since 
the field of view we consider is small (5x5 square degrees) , 
the radial comoving distance (we shall denote it by x) can be 
approximated to y. In the standard multiple- lens plane algo- 
rithm, the distance between source and observer is divided 
into N intervals separated by comoving distance Ay. The 
matter content in each interval is projected onto lens planes 
perpendicular to y-axis. The projected density contrast of 
i-th plane is defined by, 



5? IOJ (a;i,:E2) = 



dyS(xi,x 2 ,y), 



(1) 



where S = pj p — 1, yi denotes y-position of i-th lens plane. 
The two-dimensional deflection potential of i-th lens plane 
is related to the projected density contrast via the two- 
dimensional Poisson equation by 



V 2 <f l (x u x 2 ) = 



8nGp proj 



3c 2 



Sf T ° 1 (x 1 ,x 2 ) 
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Table 2. Caracteristics of PM simulation boxes. 



Model 




b 


Redshift 


range 




h _1 Mpc 


h 10 Mq 




^max 


SCDM 


80 


0.85 


0.0 


0.055545 




80 


0.85 


0.055545 


0.11585 




80 


0.85 


0.11585 


0.18147 




80 


0.85 


0.18147 


0.25306 




80 


0.85 


0.25306 


0.33136 




120 


2.9 


0.33136 


0.46332 




120 


2.9 


0.46332 


0.61591 




160 


6.8 


0.61591 


0.85950 




200 


13 


0.85950 


1.2500 




240 


23 


1.2500 


1.9054 




280 


36 


1.9054 


3.1088 


OCDM 


80 


0.25 


0.0 


0.055010 




80 


0.25 


0.055010 


0.11355 




80 


0.25 


0.11355 


0.17589 




80 


0.25 


0.17589 


0.24235 




80 


0.25 


0.24235 


0.31326 




120 


0.86 


0.31326 


0.42880 




120 


0.86 


0.42880 


0.55664 




160 


2.0 


0.55664 


0.74917 




200 


4.0 


0.74917 


1.0323 




240 


6.9 


1.0323 


1.4508 




320 


16 


1.4508 


2.1935 




440 


12 


2.1935 


3.7702 


ACDM 


80 


0.25 


0.0 


0.053997 




80 


0.25 


0.053997 


0.10942 




80 


0.25 


0.10942 


0.16641 




80 


0.25 


0.16641 


0.22515 




80 


0.25 


0.22515 


0.28581 




120 


0.86 


0.28581 


0.38084 




120 


0.86 


0.38084 


0.48137 




160 


2.0 


0.48137 


0.62541 




200 


4.0 


0.62541 


0.82485 




240 


6.9 


0.82485 


1.1001 




280 


11 


1.1001 


1.4870 




360 


23 


1.4870 


2.1341 




440 


42 


2.1341 


3.2942 



a The short side length of boxes in comoving units, the long side 
length is 2L in all cases. 
b Particle mass. 



= 3fi r 



(?) 



<5f roj (xi,x 2 ) 



(2) 



The light ray position on n-th lens plane of a ray with im- 
age positions 9 1 are computed using the multiple lens-plane 
equation, 



e n = e 1 - 



n—l 

\ ^ f{Xn 



Xi) 



*r~i a(xi)f(Xn, 



(3) 



where f(x) denotes the comoving angular diameter 
distance, defined as /(%) = K _1/ ' 2 sin _ff 1//2 x, X> 
(-K)~ 1/2 smh{-K) 1/2 x for K > 0, K = 0, K < 0, respec- 
tively, where K is the curvature which can be expressed as 
K = (Ho / c) 2 (ttm + — 1), and V± denotes either d/dxi or 
d/dx2- Note that the spatial position x 1 on i-th plane is re- 
lated to the angular position by x' = The evolution 
equation of the Jacobian matrix, which describes deforma- 



tion of an infinitesimal light ray bundle, is written by 

f(Xi)f(x 



A- 71. 



E 



a(Xi)f(Xn) 



-TJ A 



(4) 



where I denotes the identity matrix, and Ui is the, so-called, 
optical tidal matrix defined by 



(•») 



with commas denoting the differentiation with respect to x. 
It should be noticed that the first and second derivatives of 
^ in equations Q and (Q) must be evaluated at ray positions 
computed by the lens equation (|^). The Jacobian matrix is 
usually decomposed into, 



A = 



- k — 71 
-72 + w 



-72 - w 

1 — K + 71 



(6) 



where k represents convergence, I7I = (7J + 72) 1//2 is the 
amplitude of shear of light ray bundle, and u) is a net rotation 
of the beam. 

We adopt a fixed interval between lens planes by Ay = 
80/i _1 Mpc (for this choice, the long side length of simulation 
boxes are chosen so that they are multiples of 80/i~ 1 Mpc). 
The procedure to trace light rays through iV-body data 
with the multiple lens-plane algorithm can be described 
as follows: (i) computing the projected density contrast, 
^proj £ rom p ar ticle distribution in iV-body data, (ii) comput- 
ing the two-dimensional deflection potential, ^ % via Poisson 
equation (Q) and evaluating the first and second derivatives 
of it. (iii) evaluating the distance combination for each lens 
plane and performing the summation in equations ([]) and 
(^). In the rest of this section, we describe each step in some 
detail. 

(i) A simulation particle located at j/-coordinate be- 
tween yi-i and yi is projected onto the i-th lens plane at 
yt. This projection is done in parallel to the third axis, 
and thus the particle's (3:1,3:2) positions are not changed. 
Since the A^-body data is periodic in x directions^], the 
projected particle distributions are also periodic. The pro- 
jected density field is computed on a 512 2 square lattice 
from the projected particle distribution, using the triangular 
shaped cloud (TSC) assignment scheme (Hockney & East- 
wood, 1988). In this step, the most important point is the 
choice of both the smoothing scheme and the grid size. Our 
main consideration when choosing them was to maintain 
the resolution provided by the Af-body simulation and at 
the same time to remove shot noise due to discreteness in 
A^-body simulations by a relevant smoothing. We tested the 
smoothing scale by varying the size of grids, 256 2 , 512 2 and 
1024 2 . In the case of 1024 2 lattice, we found a white-noise 
contribution to the power spectrum of the lensing conver- 
gence at small scales, which is a typical signature of the 
discreteness effect (Jain et al. 2000). However, the discrete- 
ness effect was smoothed sufficiently well in both 256 2 and 
512 2 lattice cases, so the white-noise contribution was not 
detected for those cases. On the other hand, we also found 
that the small scale power was damped significantly for the 



t A-body simulation data is not periodic in y direction, because 
we used the light-cone output. 
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Table 3. Exact values of the source redshifts. In the text, they 
are refereed to as z 3 = 1, 2 and 3. 





SCDM 


OCDM 


ACDM 




1.0025 


1.0323 


1.0034 


22 


2.0422 


1.9836 


1.9721 


23 


3.1088 


2.9501 


3.0419 



256 2 case, which indicates too much smoothing. We also 
tested another smoothing scheme, cloud-in-cell (Hockney & 
Eastwood, 1988), with the same three grid sizes. In these 
cases, we found a non-negligible white-noise contribution. 
We, therefore, decided to adopt TSC with 512 2 lattice. 

(ii) Poisson equation (0) is solved to compute using 
the fast Fourier transform method with the periodic bound- 
ary condition. The first and second derivatives of VP 8 are eval- 
uated on the lattice points using the usual finite difference 
method (see Appendix A of Premadi, Martel & Matzner 
1998 for explicit expressions). 

(iii) 512 2 rays are traced backward from the observer 
point. The initial ray directions are set on 512 2 grids, which 
correspond to pixels of angular size 5°/512 ~ 0.59 arcmin. 
For each ray, we first computed ray positions on all lens 
planes in an iterative manner, using the lens equation (|^). 
The first and second derivatives of on a ray position are 
linearly interpolated from four nearest grids on which they 
were pre-computed (step (ii)). Finally, the summation in 
equation (^) is performed. The light ray positions and four 
components of the Jacobian matrix on desired source planes 
are stored. Exact values of the source redshifts used in this 
paper are summarized in Table [j. 

We performed 40 realizations of the underlying density 
field for each cosmological model by random shifts of the 
simulation boxes in the {x\,X2} directions using periodic 
boundary condition. Note that lens planes coming from the 
same box are shifted in the same way to maintain the clus- 
tering of matter in the box. It is important to note that 
the 40 realizations are not rigorously independent because 
they come from the same set of the simulation boxes. It is 
expected that the cosmic variance can not be correctly in- 
cluded because of this lack of independence, but this does 
not change our results because on angular scales of our in- 
terest (below 100 arcmin) it is very small. 



3 NUMERICAL RESOLUTION 

The angular resolution of the ray-tracing simulation is ba- 
sically limited by the spatial resolution of the TV-body sim- 
ulations. Since we used PM iV-body code, the spatial res- 
olution is simply r rcs ~ L/256. This can be converted into 
the angular scale using the angular diameter distance rela- 
tion, r rcs = f(x)0 ICB . The bottom panel of Figure]^ shows 
such angular scales computed for three cosmological mod- 
els. Thanks to the tiling technique, 9 rcs is almost constant 
(dies ~ 1.5) at redshifts higher than 0.2. It increases as 
#rcs oc x" 1 at z ^ 0.2, because of the constraint on the 
minimum box size. The impact of this worse resolution on a 
measurement depends on the quantity one wants to exam- 
ine. Roughly speaking, quantities which come from solving 




12 3 

z 

Figure 2. Top panel: Distance combination appearing in the 
equation (0) for z s = 1 and 3 as a function of lens redshift. Middle 
panel: Distance combination appearing in the equation (Q) (nor- 
malized by Hubble length, c/Hq) for z s = 1 and 3 as a function 
of lens redshift. Bottom panel: 8 res as a function of redshift (see 
text) . 



the lens equation (|3|) are more influenced by the worse res- 
olution than ones related to the Jacobian matrix obtained 
via equation (Q). 

The lensing excursion angle plays a crucial role on stud- 
ies of the weak lensing effects on the CMB and is denned 
by the difference in deflection angles between nearby two 
light rays. Each light ray trajectory is computed by the lens 
equation. As the expression, eq. (0), indicated, the distance 
combination f(Xn — X<)//(X») a °t s as a weight function. 
The top panel of Figure ti shows this distance combination 
as a function of the lens redshift for two cases, z s = 1 and 
3, and clearly indicates that lensing deflections due to struc- 
tures at low redshifts are more weighted than those by high 
redshift structures. It can be, therefore, expected that the 
worse angular resolution at low redshifts has a significant 
impact on the effective accuracy of the light ray trajectory. 
In order to check it, we compute the variance in the lensing 
excursion angles of the two nearby light rays which would be 
observed with a separation 9 if there was no lensing. The de- 
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1 10 100 1 10 100 1 10 100 

0[arcmin] 0[arcmin] 0[arcmin] 

Figure 3. Top: The root-mean-square of lensing excursion angles 
between nearby light ray pairs divided by their intrinsic separa- 
tion 0. 5a is either Scti or <5»2j and the averaged value is plotted. 
The open circles are for measurements of z s = 1, while the filled 
triangles are for z s = 3. The error bars denote the standard devi- 
ation computed among 40 realization. The solid and dotted lines 
represent the nonlinear and linear prediction, respectively. Bot- 
tom: Ratio of the measurements to the nonlinear prediction for 
Zs = 3. 



flection angle of a light ray is simply the difference between 
its angular position in the first lens plane and that in the 
source plane, i.e., a = 9 1 — 6 n . The lensing excursion angle 
between nearby rays A and B is simply 3a — a A — a B , 
and their intrinsic separation is — \0 n,A — 6 n ' B \. In Fig- 
ure^, the results are plotted together with the semi-analytic 
predictions by the power spectrum approach (Seljak 1994; 
1996). The measurements agrees very well with the nonlin- 
ear semi-analytic prediction at scales larger than 10 arcmin. 
Below that scale, the measurements are depressed reflecting 
the effective resolution limit. We may, therefore, conclude 
that the angular range, where our analysis related to the 
lensing excursion angle is reliable, is between 10 arcmin and 
2 degree. 

Although we are not dealing with statistics of the lens- 
ing convergence in this paper, it is instructive to compare 
the effective resolution of the lensing convergence with that 
of the lensing excursion angle. The lensing convergence is 
computed by solving the evolution equation of the Jacobian 
matrix, eq. (^) . The distance combination appearing in this 
equation is f{Xi)f(Xn ~Xi)/ '/(x«)> an d is shown in the mid- 
dle panel in Figure fel This distance combination has a peak 
at intermediate redshift depending on the source redshift. 
Therefore, the worse resolution at low redshift lens planes 
has only small impact on the effective resolution of lens- 
ing signals, k, and 7, except for cases of very low source 
redshifts (z s 0.5). Figure ^ shows that the two-point 
correlation function of the lensing convergence defined by 
C K (9) — {/t(</>)«(0 + 6)) compared with the semi-analytic 
prediction (Bartelmann & Schneider 2001). The measure- 
ments agree well with the predictions down to about 2 ar- 
cmin. Therefore, the effective resolution of the lensing con- 
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1 10 100 1 10 100 1 10 100 

0[arcmin] 0[arcmin] 0[arcmin] 

Figure 4. Same as Figure^ but for two-point correlation function 
of the lensing convergence. 




Figure 5. The probability distribution function of the lensing 
excursion angle normalized by its intrinsic separation (the solid 
curves). Top panels are for z s = 1 and bottom panels for z s = 3. 
Left and right panels are for light ray pairs with intrinsic sepa- 
ration 8' < 6 < 12' and 20' < 6 < 32', respectively. The dotted 
curves show Gaussian distribution with the a computed the mea- 
sured PDFs 



vergence is found to be much better than that of the lensing 
excursion angle. 
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4 STATISTICAL PROPERTIES OF LENSING 
LENSING EXCURSION ANGLES 

Figure |E] shows the probability distribution function (PDF) 
of the lensing excursion angles normalized by its intrinsic 
separation. Since the vector field 8a{<p) has no special di- 
rection, we take both two components, Sai and Sa2, to com- 
pute the PDFs. ft is clearly shown in Figure ^ that the 
PDFs consist of two contributions: A Gaussian distribution 
at inner part, and the exponential tail at the outer part. 
The inner part of the PDFs are fitted reasonably well by a 
Gaussian distribution with the a computed from the PDFs 
(root-mean-square computed from the PDFs was used for 
a). We should here noticed that the measured standard de- 
viation agrees very well with the semi-analytical prediction 
as shown in Figure |S| 

The fact that the PDFs consist of two contributions 
suggests that there are two different processes that make the 
lensing excursion angle. One, which makes Gaussian distri- 
bution, might be secular small (random) deflections due to 
linear density fluctuations along each light ray path. This 
is explained by the fact that the lensing deflection angles, 
which are due to Gaussian random fluctuations (such as the 
linearly evolved density fluctuation field) , are Gaussian ran- 
dom field. As the light ray travels longer distance, the ray 
can undergo more fluctuations. Therefore the width of Gaus- 
sian distribution increases as the source redshift becomes 
higher. The other, which makes the exponential tail, might 
be a single (or possibly multiple) coherent scatter by a non- 
linear structure such as a galaxy or a cluster of galaxies^} 
As the separation of a light ray pair decreases, smaller scale, 
strongly nonlinear structures can contribute to a coherent 
scatter. Therefore the exponential tail becomes more promi- 
nent for a small separation case than a larger one as shown 
in Figure |[ 

Does the exponential tail make a significant influence 
on studies of weak lensing effects on the CMB ? As far as 
the weak lensing effects on the CMB power spectrum con- 
cerned, it has no effect, because the crucial assumption that 
the lensing excursion angle is (in a statistical sense) much 
smaller than the intrinsic separation is true even in the pres- 
ence of the exponential tail. This can not be directly demon- 
strated by our simulations as the light rays are not followed 
up to LSS, but can be proved by the semi-analytic predic- 
tion which tells that the root-mean-square of the lensing 
excursion angle is much smaller than the intrinsic separa- 
tion angle. The validity of the semi-analytic prediction was 
supported by our numerical simulation (Figure |^). 

Figure n shows the Kurtosis of the PDF of 8a defined 
by (See 4 ) / (5a 2 ) 2 — 3. The Kurtosis decreases as the source 
redshift becomes higher. On scales larger than fO arcmin 
(below that, the lensing effects on the CMB will be hardly 
detected) the Kurtosis less than unity at z 3 = 3, and thus, 
at the LSS, it must be smaller than that. It can be, there- 
fore, said that the exponential tail or, more generally, the 



S Computations of a lensing optical depth of strong lensing events 
by a galaxy or a clusters of galaxies as well as a small probability 
of finding multiple imaged QSOs in a QSO catalog tell us that 
light rays which undergo strong lensing more than twice times 
are very rare (e.g., Chiba & Futamase 1999; Hattori, Kneib & 
Makino 1999, Hamana, Martel & Futamase 2000). 
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Figure 6. Kurtosis of lensing excursion angle 8a, that is, 
(8a 4 ) / {8a 2 ) 2 — 3, as a function of the intrinsic separation of light 
ray pairs. Filled circles, filled squares and filled triangles are for 
Z s = 1, z s = 2 and z s = 3, respectively. 



non- Gaussian nature in the PDF is very unlikely to have a 
significant effect on the weak lensing of the CMB. 



5 SUMMARY 

We performed ray-tracing simulation combined with sets of 
large iV-body simulations. The use of the tiling technique 
enable us to explore a wide angular dynamic range with the 
efficient PM iV-body code. The angular resolution is lim- 
ited by the spatial resolution that A-body simulations have. 
Since the minimum box size is constrained to keep large scale 
modes that may contribute to the weak lensing, the angu- 
lar resolution of the low redshift boxes becomes worse. We 
found that the tiling technique particularly suits for study 
of the cosmic shear statistics, because the worse angular 
resolution does not make a serious impact on an effective 
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angular resolution of the lensing convergence and shear. On 
the other hand, this worse resolution make a stronger effect 
on the effective resolution of the lensing excursion angles. 
One possible solution to improve the resolution is to use a 
higher resolution simulation such like P 3 M code for a few 
lowest redshift boxes. 

We have numerically examined statistical properties of 
the lensing excursion angles. We found that the variance in 
the lensing excursion angles predicted by the power spec- 
trum approach is in good agreement with our numerical 
results. We found that the PDF of the lensing excursion 
angles is not simply Gaussian but has an exponential tail. 
However, it can be safely concluded that the exponential 
tail, or more generally non-Gaussian nature in the PDF has 
no significant effects on the weak lensing of the CMB. These 
results demonstrate a validity of using the power spectrum 
approach to compute lensing effects on the CMB. Further- 
more our results support the validity of a simple numerical 
simulation method to obtain a lensed CMB map in which 
both the lensing displace filed and the CMB temperature 
map are generated assuming the Gaussian statistics and the 
lensed CMB map is obtained by a mapping (e.g., Takada & 
Futamase 2001). 
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